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1.0 SUMMARY 


This document describes the design and usage of a pilot program for calculating pressure 
distributions over harmonically oscillating airfoils in transonic flow. The procedure is based 
on transonic small perturbation theory, in which the velocity potential is separated into 
steady and unsteady parts, and the resulting differential equation for the unsteady flow is 
linearized by assuming small disturbances. The steady velocity potential distribution, which 
must be obtained from some other program, is required for input. Time is not a variable in 
the program since harmonic motion is assumed. The differential equation for the unsteady 
complex velocity potential is linear, with spatially varying coefficients that are dependent 
on the steady velocity potential distribution. The numerical solution is obtained through a 
finite difference formulation and either a line relaxation solution or a full direct solution. 



2.0 INTRODUCTION 


This document describes the design and usage of the FORTRAN IV digital computer 
program A344. This pilot program was written as part of a continuing research effort to 
develop a method for analyzing the transonic flow about harmonically oscillating airfoils 
and wings. 

The program uses a procedure based on small perturbation theory in which the velocity 
potential is separated into steady and unsteady parts. The resulting equation for the 
unsteady velocity potential is linearized by assuming small perturbations. The solution is 
obtained using finite difference techniques. Boundary conditions are applied on the outer 
mesh boundaries, the plane of projection of the airfoil (i.e. on the slit y=0), and over the 
wake (i.e. the plane vortex sheet downstream of the trailing edge). The derivation of the 
equations and the application of the finite difference methods are discussed in detail in 
references 1 through 5. The program computes the unsteady velocity potential and the 
resulting unsteady pressure distributions. It requires as input the potential distribution from 
a steady state transonic small perturbation program (e.g. ref. 6). Conservative differencing 
is used for subsonic points and nonconservative differencing for supersonic points and across 
shocks. The program uses a rectangular grid. 

The solution may be obtained by using either a relaxation procedure or a direct procedure. 
The procedure and routines for the direct solution, which is out-of-core, is presented by Yip 
in reference 7. 

Features of this program include options to: 

• Analyze sections with or without thickness 

• Use row or column relaxation, or use a direct solution 

• Take advantage of antisymmetry for symmetric airfoils at zero angle-of-attack to 
perform the solution over half the finite difference space 

• Analyze flows with either subsonic or supersonic freestreams 

• Obtain the pressure distributions for vertical translation, pitch and control surface 
motions 

A parallel program for three-dimensional flow but using only relaxation procedures is dis- 
cussed in reference 8. 
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3.0 SYMBOLS 


b Semichord 

I Column index 

lA Index of first column aft of control surface hingeline 

IMAX Number of columns in mesh grid 
10 Index of first column aft of airfoil leading edge 

II Index of column at trailing edge or immediately aft of trailing edge 

J Row index 

JM Index of row immediately below airfoil 

JMAX Number of rows in mesh grid 
K Transonic parameter, (l-M^)/(M^e) 

M Freestream Mach number ■ , 

PHIO Steady-state velocity potential array 

PHI 1 Unsteady velocity potential array 

U =K-(7+ IXd^^/dx) 

X,Y Scaled coordinates, X = Xp, Y = pYp 

Xp, Yp Physical coordinates, made dimensionless with semichord, b. 

7 Ratio of specific heats 

5 Thickness ratio or measure of camber and angle of attack 

e (5/M)2/3 

Xj coM/(1-m2) 

p Scale factor on Yp, p = 

'Pq Steady perturbation velocity potential 

Unsteady perturbation velocity potential 

CO Angular reduced frequency (semichord times frequency in radians per second 

divided by freestream velocity) 
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4.0 DISCUSSION 


The equations for the relaxation portion of the program of this document are fully derived 
and presented in the Appendix of reference 1 . The notation used in the program is for the 
most part taken directly from that report. The geometry used for the program is shown in 
figure 1 . The direct solution portion of the program differs mainly in the form of the outer 
boundry condition where out-going wave boundary conditions are used rather than Klunker- 
type boundary conditions. 

The program includes four solution options: 

1 Relaxation, full-space, subsonic freestream 

2 Direct solution, half-space, subsonic freestream 

3 Direct solution, half-space, supersonic freestream 

4 Direct solution, full-space, subsonic freestream 

The relaxation procedure is an iterative process in which the solution is obtained as sequence 
of relatively small problems. In the direct solution, the difference equations are solved all at 
once. Generally, relaxation solutions use significantly less core storage and more central 
processing time. The direct solutions use more core storage and less central processing time. 
However, as noted in references 2 and 3, relaxation solutions diverge for values of above 
a certain critical value. This divergence severely limits the useful range of relaxation pro- 
cedures, particularly in terms of reduced frequency at higher values of Mach number. The 
estimation of the critical value of Xj , which is a function of Mach number, reduced 
frequency, and the dimensions of the solution space, is discussed in section 4.0 of reference 2. 

The equations of fluid flow (equations 24 and 26 of ref. 1) are set up for row relaxation as: 
SUB1(I,J) * <^i(I-2,J) -t- SUB (I,J) *(pi(I-l,J) + DIAG(I,J) * v?ld,J) 

-I- SUPER(I,J) *(pi(I+l,J) = RHS(I,J) 

The right-hand side term, RHS(I,J), contains values from both the current and previous 
iterations. For example, for row relaxation proceding from lower to upper boundary, 
RHS(I,J) includes values from the current iteration for the row below the Jth row (row J-1) 
and values from the previous iteration for the row above the Jth row (row J+1). For 
column relaxation, the unknowns on the left-hand side of the equations are 
on the left-hand side of the equations are ip]^(I,J-l), sPjfl,!) and i^j(I,J-i-l). For the direct 
solution all five are on the left-hand side. 

Solutions for the direct procedure are obtained using an out-of-core solution routine. The 
procedure and corresponding routines are presented in detail by Yip in reference 7-. 

The program calculates the modal input data for vertical translation, pitch and control 
surface motions. A separate fun is required for each mode. 
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(Notation matches that of NASACR-2257) 


Figure 1 — Solution Mesh for Two-Dimensional Problem 
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5.0 COMPUTER PROGRAM USAGE 


5.1 MACHINE REQUIREMENTS 

A344 executes on a CDC CYBER 175 or similarly compatible computer. 

5.2 OPERATING SYSTEM 

A344 was designed for a KRONOS 2.1 or NOS 1.2 operating system. 

5.3 STORAGE ALLOCATION 

Relaxation option, NVER=1 ; requires cm=l 20000 (octal words) 

Direct solution options: 

Half-space, subsonic freestream, NVER=2; cm=220000 
Half-space, supersonic freestream, NVER=3; cm=2 10000 
Full-space, subsonic freestream, NVER=4; cm=3 10000 

5.4 TIMING 

Timing is hardware- and operating system- dependent, as well as being a function of the 
number of mesh points used in the finite difference grid. 

The following direct solution examples were run on a CDC Cyber 1 75 with NOS 1 .2 operat- 
ing system. 

1 . Full-space solution, 56 x 60 mesh; 70 CPU Sec. 

2. Half-space solution, 64 x 46 mesh; 1 2 CPU Sec. (a 64 x 23 mesh is actually calculated) 

5.5 FILE I/O 

The program card of A344 is as follows: 

Program TEA344 (INPUT= 102, OUTPUT=1002, TAPE7, TAPE5=INPUT, 
TAPE6=OUTPUT,TAPEl =1 002,TAPE1 0=1 002,TAPE98=5 1 2, 

TAPES) 

The buffer sizes are given in decimal form, with the default value being 1027 decimal 
(2003g). 

5.5.1 FILE UTILIZATION 

The file named TAPEl contains one steady state velocity potential distribution (PHIO) 
from a separate program. This matrix must be present for a wing with thickness (see ref. 3 
and input variable MSTST). 
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The file named TAPEIO is the file where the PHIl matrix is stored. If the relaxation solu- 
tion procedure is used, the starting PHIl (if available) is stored on TAPEIO. After a success 
ful execution of the program (using either the relaxation or direct solution procedures), 
TAPEIO will contain the new PHIl matrix. If an old PHIl was stored on TAPEIO, it is 
replaced by the new PHIl. The format for TAPEIO is described in section 5.5.2. 

5.5.2 FILE FORMATS 

Binary files, TAPEl and TAPEIO. 

The files TAPEl and TAPEIO are binary I/O files. The data on TAPEl is real and includes 
the number of rows and columns, the x and y coordinates, and the steady state velocity 
potential. TAPEl is read with the following sequence of READ statements: 

READ(l) IMAX,JMAX 

READ (1) (X(I), 1=1, IMAX), (Y(J), J=l, JMAX) 

READ(l) ((PHI0(I,J),I=1,IMAX), J=l, JMAX) 

The PHIO matrix is required only if the airfoil has finite thickness, in which case set 

MSTST = 0. The parameters IMAX and JMAX, and the coordinate arrays X and Y are 
required in all cases. 

The data on TAPEIO is complex with a real- and imaginary-part for each element of the 
unsteady velocity potential matrix, PHIl. TAPEIO is read with the statement 

READ(l) ((PHI(I,J),I=1,IMAX), J=l, JMAX) 

If the relaxation solution is being used, a previously calculated PHIl may be used as a 
starting point for the new calculations by setting INC^O. 

As mentioned previously the complex binary files will contain a real part and an imaginary 
part for each element of the PHIl matrix. 

BCD Files 

The BCD files INPUT, OUTPUT, and TAPE? follow the standard FORTRAN and system 
formats for that type of file. 

BCD files are those which deal with character printing or reading. A344 has three of these; 
INPUT, OUTPUT, and TAPE?. 

INPUT, also called TAPE5, is the file which contains cards or card images. Program card 
(or card image) input is fully described in section 5.7. 

OUTPUT, also called TAPE6, is the file on which the program places the primary printed 
information, (see sec. 3.8 and appendix) 


7 



TAPE? is also a print file. The user may disregard it unless he is executing the program on 
a terminal where the primary OUTPUT print file is to be printed later. It will print a 
summary during execution telling the user how convergence is proceeding. 

Usage of TAPE? is EGO, IN, OUT, OUTPUT, (see sec. 5.6) 

TAPE? was primarily used in development of the program. Terminal usage of the program 
should be limited as terminal execution is usually very expensive. 

5.6 CONTROL CARDS 

The following control cards can be used to load and execute TEA344 from permanent file: ■ 

JOBN, T20, CMS 10000, POO. 

USER, ACCTNO, PASWD. URNAME/PN/MS/ORG 
GET, XXPROG=TEA344. 

GET, TAPE 1=URPHI0 
GET. XXDATA-URINPT. 

XXPROG (XXDATA). 

If the relaxation solution version is to be used (NVER=1) and an existing PHIl distribution 
is to be used as a starting point, a GET, TAPE10=URPHI10 card should be added, if PHIl 
is to be saved, a SAVE, TAPEIO = NEWTPIO control card is needed. 

5.1 PROGRAM INPUT 

5.?. 1 GENERAL REMARKS 

The input to A344 is of two forms, disk file/tape input (binary input) and card input 
(BCD input). Disk file/tape input may consist of input point locations and the PHIO and/or 
PHIl distribution. An input PHIl distribution is indicated if a previously calculated PHIl 
matrix is to be used to start the iteration process. If the user is starting an iteration pro- 
cedure from scratch, there is no PHIl input and the initial PHIl distribution is all zeros. 

The subroutine INCOND generates these zeros if MSTS = 0. There is also no PHIl input if 
a direct solution procedure is used. PHIO is the steady-state distribution from another pro- 
gram. If a flat plate solution is sought, PHIO would not be input. 

A description of disk file/tape formats is given in section 5.5.2 and a listing of the input for 
a sample problem is presented in the appendix. 

The card input consists of field dependent input and namelist free field input. The field- 
dependent input is defined in the format column of table 1. as a spceific field (i.e., FI 0.2, 
AlO, 15). The namelist data will be represented in the same column by the namelist name 
“PARAM.” 
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It should be noted that the namelist input set is not the same for the four solution options. 
The input variables ILAX, CONPXT, and CONE6 are read only for the relaxation solution, 
version 1. The input variables LW and MXRR are read only for the direct solutions, versions 
2, 3 and 4. Input variably ICKPRI is read only for version 2. 

Some of the features of namelist input are: 

1 . Card field consists of columns 2 through 80. 

2. List consists of a $ list-name in column 2, in our case, $ P ARAM, followed by a series of 
specifications continued on as many cards as required and terminated by a $ END. 

3. Specifications are of the form: 

a. Vname = Value 

b. Where Vname is an array, Vname = Valuel ,...,Valuen. Where Vname is one of 
the variable names for the list, value is the associated value(s). Value may be an integer, 
a floating point number in normal or exponential form, or in the case of a logical 
variable (specifically the options), of the form. 

- .T. or .True . indicating true or on 

- .F. or .False, indicating false or off 

4. Specifications must be separated by commas. There is no comma between the last 
specification and terminating $. 

5. Embedded blanks are allowed except within the $ PARAM , variable name, or value. 

At least one blank must separate the $ list-name and the first specification. 

6. The order of appearance of variables on the card(s) is not important; the spelling is. 

7. Any or all of the variables may be left out of the list, e.g., $ list-name, $ is legitimate. 
This assumes, of course, that there is a legal default value associated with the variable(s) 
not included in the list. 

5.7.2 LIMITATIONS 

The following are size limitations within the program. 

• 3<IMAX<75 IMAX = Maximum number of X - nodes (in flow direction). 

• 3 < JMAX < 60 JMAX = Maximum number of y - nodes (in vertical direction). 

• 0 < OMEGA 

• FSMACH < 1 For subsonic freestream option. 

Note: The PHIO and PHIl distributions must also correspond to the limitations on 

the XY mesh. 
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5.7.3 DATA STACKING 

Note: All coordinates and lengths, except for the location of the pitch axis, 

are entered as scaled quantities, 

X = X /b 
Y=M*Yp/b 

where b is the semi chord and the subscript p means physical coordinate. The pitch axis is 
entered as Xp / (2b) units aft of leading edge. 

The scale factor p, as well as the parameter e, is undefined for the flat plate for which 5=0. 
The program is coded such that it will “blow up” if 5 = 0 is used as input. Therefore, a 
finite value of 5 should be used for flat plate configurations. The unsteady velocity potential 
distribution will be scaled accordingly; however, the pressure coefficient is independent of 
6 since it is found by multiplying the x - derivative of the potential by the factor e. 
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Card number 
1 


2 

3to N 


Table l.—User Input Variables 

Variable name Format Description 

VERSION 7X,I2 Denotes version of solution to be used. 

VERSI0N=1 uses the full-space relaxation. 
VERSI0N=2 uses haif-space, direct solution, subsonic 
freestream. 

VERSI0N=3 uses half-space, direct solution, 
supersonic freestream. 

VERSION=4 uses full-space, direct solution, 
subsonic freestream. 


Title 

8A10 

80-character title for this problem 

FSMACH 

PA RAM 

Mach number 

DELTA 

PARAM 

Thickness ratio 

THETA 

PARAM 

Amplitude of the oscillatory angle of attack in degrees 

OMEGA 

PARAM 

Angular reduced frequency 

GAMMA 

PARAM 

Ratio of specific heats for flow medium 

AL 

PARAM 

Chord of control surface 

IMAX 

PARAM 

Maximum X node count in users mesh 

JMAX 

PARAM 

Maximum Y node count in users mesh 

NVOL 

PARAM 

No. of far-field updates before the volume integration is 
included in the far-field calculation. Limited experience 
indicates that the volume integration does not appear signif- 
icant to results. The volume integration may also be sup- 
pressed by setting IS = 0. 

IS 

PARAM 

Starting X node limit for volume integral for the wing; 
currently not used. Intended only for relaxation 
solutions, if IS=0, the volume integral is not included in the 
farfield evaluation. 

IE 

PARAM 

Ending X node limit for volume integral for the wing; 
currently not used. Intended only for relaxation solutions. 

JS 

PARAM 

Starting Y node limit for volume integral for the wing; 
currently not used. Intended only for relaxation solutions. 

JE 


Ending Y node limit for volume integral for the wing; 
currently not used. Required only for relaxation solutions. 

NMAX 

PARAM 

Maximum number of iterations to be allowed without 
convergence. Required only for relaxation solutions. 

NA 

PARAM 

Far-field update cycle control; updates the far-field each NA 
iteration. Required only for relaxation solutions. 

ERROR 

PARAM 

Error difference. When the maximum difference between 


PHI1 distributions of consecutive iterations is less than 
ERROR, the program stops iterating. Required only for 
relaxation solutions. 
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Table 1.— (Continued) 


Card number Variable 

Format 

Description 

NP 

PARAM 

Prints pressure distribution every NP iterations. Required only 
for relaxation solutions. 

INC 

PARAM 

Restart variable. IF INC^0,start with the PHI1 distribution 
on TAPE10; if INC=0, start with a PHI1 distribution of zeros. 
Required only for relaxation solutions. 

ORF 

PARAM 

Overrelaxation factor used for subsonic nodes to accelerate 
convergence,1^RF<2. Required only for relaxation 
solutions. 

URF 

PARAM 

Underrelaxation factor used for supersonic nodes to accelerate 
convergence, 0 <URF<1. Required only for relaxation 
solutions. 

MSTST 

PARAM 

When MSTST=0, a PHIO distribution from a steady-state 
solution is read from TAPE1 and used to calculate U(I,J). 

When MSTST 7^0, U(1,J) is set to K and thickness effects are not 
included in the analysis. This is the f/af p/afe analysis. 

ISWEEP 

PARAM 

This variable along with ILAX determines the sequence in 


which line relaxation solutions are obtained. With I LAX = 0 
(row relaxation): (1) if ISWEEP = 0, rows will be solved from 
the upper and lower boundaries, alternating in toward the 
wing; (2) if ISWEEP = 1 , rows will be solved from the airfoil, 
alternating out to the upper and lower boundaries; (3) if 
ISWEEP = 2, rows will be solved consecutively from lower 
boundary to upper boundary. With I LAX = 1 (column 
relaxation): (1) if ISWEEP = 0, columns will be solved start- 
ing at the trailing edge of the airfoil and moving forward to 
to the upstream boundary, then going back to the trailing edge 
and moving aft to the downstream boundary; (2) if ISWEEP 
= 1 , columns are solved from upstream boundary to down- 
stream boundary. Required only for relaxation solutions. 
Default values are I LAX=0 and 1SWEEP=0 

I LAX PARAM With I LAX = 0, the line relaxation solution procedure uses 

rows of points, the points forming a iine parallel to the airfoil. 
With I LAX = 1, the line relaxation procedure works with 
columns of points which are perpendicular to the airfoil. 
Required and read for the relaxation solution only. 

CONPXT PARAM Constants required for convergence of row relaxation (refs. 

2 and 3) Read only when NVER = 1. 

Ip Not used. 

MXRR PARAM Maximum number block rows. Required for direct solutions 

only. Default value Is 29. MXRR must not be greater than 29. 
Not read when NVER = 1 . 
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Card number Variable 
ICKPR1 PARAM 


LW PARAM 


MODE PARAM 

PAXIS PARAM 


Table 1 .—(Concluded) 

Format Description 

Checkprint option parameter available only for half-space 
subsonic freestream version (NVER=2). If ICKPRIj^O, 
a certain limited amount of the intermediate calcula- 
tions is printed out for check purposes. This includes the 
NDBLK array, the IN matrix and several of the block 
coefficient matrices, WW. If ICKPR1 = 0, this data is not 
printed. 

This is the LWORK parameter required by the out-of-core 
direct solution procedure, (see ref. 7). The default value is 
matched to the maximum permissible values of I MAX and 
JMAX. (see section 5.7.2) 

Selects mode shape. Use MODE = 1 for vertical translation, 
MODE = 2 for pitch, and MODE = 3 for control surface 
motion. 

Location of pitch axis measured aft from leading edge. 

Note that units are Xp/(2b) 
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5.8 PROGRAM INPUTS 


5.8.1 PROGRAM RESULTS 

A listing of the output for a sample problem using the direct solution is presented in the 
Appendix. 

The printed output of the program consists of an initial block of information printing back 
the user’s input followed by information identifying some of the program options the user 
selected. 

The X and Y mesh data come next, as read from either cards or binary file (TAPE!). The 
mesh data are followed by the values of calculated variables, and the FU array which is the 
downwash on the upper surface of the wing. 

The parameter ICKPRl has been set to 1, resulting in the printing of the NDBLK array 
(the sizes of the matrix blocks), the IN matrix (the blocks with nonzero elements are 
identified with integers), values of MXR and MXC, and several coefficient block matrices. 

Next is printed a summary of the augmented matrix, and an indication of the quality of the 
solution as measured by a “norm” as defined in reference 6 and the largest residual. For 
example, “WW BLOCK 1=4” h the block matrix resulting from putting the block matrices 
6,7,8 (see IN matrix) together, i.e. [6|7|8] . 

Finally, the jump in pressure coefficient across the wing is printed. 

A second sample problem using a relaxation solution rather than a direct solution is also 
included in Appendix A. The initial parameters and mesh data are listed first as for the 
direct solution. Next the ERROR is printed out after each iteration. The program also prints 
out intermediate (and unextrapolated) pressure coefficients every NA iterations. The 
final interpolated pressure coefficients are printed out after the solution converges or after 
NMAX iterations, whichever occurs first. 

5.8.2 PROGRAM DIAGNOSTICS 

“SOLUTION FAILED TO CONVERGE IN ITERATIONS lERR, JERR ERRMAXl 

^ ^ ” indicates the largest error found at the indicated XY node location is 

still larger than the user-specified standard and that the maximum number of iterations has 
been attained. This is found in the relaxation solution option. 

“SOLUTION CONVERGED. MAXIMUM ERROR IS ” indicates the user error stand- 

ard has been reached. This is found only in the relaxation solution option. 
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6.0 COMPUTER PROGRAM DESCRIPTION 


6.1 OVERLAY STRUCTURE 

The TEA344 program consists of a (0,0) level overlay, two primary overlays and three 
secondary overlays. 

t 

The (0,0) overlay selects whether the solution is to be a relaxation or direct procedure. 

The primary level (1,0) overlay is the relaxation solution procedure. The (2,0) overlay selects 
the particular form of the direct solution. 

The secondary overlays (2,1), (2,2), and (2,3) perform direct solutions for subsonic free- 
stream, half-space; supersonic freestream, half-space; and subsonic freestream, full-space 
conditions. 


6.2 COMMON BLOCK USAGE 

Table 2. describes usage of blocks of common variables. 
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TABLE 2.*^ Blocks of Common Variables 


Block 

BO 

CONE6 

CONST 

EPS 

FUL 

FXY 

GH 

I MESH 

IRLAXF 

IRELAX 

IVOL 

L 

MATRIX 

PHI 

SOBER 

SUBSUB 

SUM 

U 

V 

WAK 

WAKE 

XY 


Usage 

Wake integral quantities 

Row relaxation constant 

Control and calculation constants 

Mach number, angle of attack and amplitude of oscillation 

Modal data for boundary conditions on camber line 

Functions of X and Y of the flow field mesh 

Bessel function values 

X and Y mesh constants 

Array to determine whether column is all subsonic or contains supersonic points 
Integers for row relaxation options 

Volume integral constants for Klunker-type boundary conditions 
Far-field integral quantities, Klunker-type boundary conditions only 
Coefficient arrays required to set up solution matrices 
Unsteady velocity potential array, PH II 
Additional coefficient arrays required for solution matrices 
Additional coefficient array for row relaxation solution 
Far-field integral quantities, Klunker-type boundary conditions only 
Steady flow quantity, U = K-(7-H 

Terms involved in evaluating the PHI1 on outer boundary; used in relaxation 
procedure only 

Wake terms for direct solution 

Wake integral quantities for Klunker-type boundary conditions only 
X and Y locations of mesh points 


6.3 PROGRAM/SUBROUTINE DESCRIPTION 

Short descriptions for the various programs and subroutines of A344 are given in this 
section. Routines are ordered alphabetically. Maps of the overlay structures and the sequence 
in which the programs and subroutines are used are shown in figures 2 and 3. 

1 . SUBROUTINE AIRFOIL 
Calculates boundary conditions over airfoil. 

2. PROGRAM A344VB1 

Performs direct solutions for half-space, subsonic freestream condition. 

3. PROGRAM A344VB2 

Performs direct solution for half-space supersonic freestream condition. 

4. PROGRAM A344VB3 

Performs direct solution for full-space subsonic freestream condition. 

5. SUBROUTINE BESNIS 

Required for calculation of Bessel functions. 

6. SUBROUTINE BESNKS 

Required for calculation of Bessel functions. 

7. SUBROUTINE CPR 

Calculates the difference in pressure coefficients across the airfoil. 

8. SUBROUTINE DATE 
Obtain date; system subroutine. 

9. SUBROUTINE DELPHI 

Calculates the jump in velocity potential across the airfoil. The velocity potential is found at 
the wing surface by extrapolation using values at the two points below (or above) the wing 
before the jump is calculated. 
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Figure 2.— Overlay and Program Map 
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Figure 3.— Overlay, Program and Subroutine Map 



10. SUBROUTINE ECTCSM 


Subroutines from out-of-core solution systems. See reference 7. 

11. SUBROUTINE ETRDSOL 

Subroutines from out-of-core solution systems. See reference 7. 

12. SUBROUTINE ETSMGEN 

Subroutines from out-of-core solution systems. See reference 7. 

13. COMPLEX FUNCTION F 

Evaluate far field function at specified outer boundary point. 

14. SUBROUTINE FARFLD 

Calculates Klunker-type boundary conditions on outer mesh boundaries. 

15. SUBROUTINE FLWTEST 

Checks column of mesh points for supersonic flow. 

16. SUBROUTINE HANTRM 

Evaluates 3 expressions involving Hankel functions for SUBROUTINE FARFLD. 

17. SUBROUTINE HANKEL 

Evaluates the complex-valued Hankel function of first or second kind for real arguement 
X and integer order. 

1 8 . SUBROUTINE INCOND (INC) 

Provides an initial distribution of unsteady velocity potential. 

19. SUBROUTINE MATRXCF 

Calculates coefficients for column relaxation solutions. 

20. SUBROUTINE MATRXFL 

Routine to calculate coefficients for full-space, subsonic freestream solutions. 
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21. SUBROUTINE MATRXC2 

Routine to calculate coefficients for half-space subsonic freestream solutions. 

22. SUBROUTINE MATRXC3 

Routine to calculate coefficients for half-space supersonic freestream solutions. 

23. SUBROUTINE MESH 

Routine to calculate mesh functions for relaxation solutions including Klunker-type 
boundary conditions. 

24. SUBROUTINE MESH2 

Routine to calculate mesh functions for direct solutions including out-going wave boundary 
conditions. 

25. SUBROUTINE SECOND 

Determines central processor time; system subroutine. 

26. SUBROUTINE STDYST (MSTST) 

Inputs steady velocity potential and evaluates U(I,J). 

27. PROGRAM TEA344 

Selects either relaxation or direct solution. 

28. PROGRAM TEA344A 
Performs relaxation procedure. 

29. PROGRAM TEA344B 

Selects one of three direct solution procedures. 
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APPENDIX 


Input for the first sample problem is for a direct solution of an airfoil in a Freon flow. 


Card Columns 

No. 1 11 21 31 41 51 ^ 

1 VERSION 2 

2 SAMPLE PROBLEM, HALF SPACE, SUBSONIC FLOW 

3 $PARAM FSMACH=0.9, DELTA=0.06, THETA=1.5, OMEGA=0.30, 

4 GAMMA=1.135, 

5 MSTST=0, 

6 MODE=2, 

7 PAXIS=0.0, 

8 ICKPRI=1, 

9 SEND 

Output for the first sample problem starts on the next page. 
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•53597E>02 

3« 

0. 

C. 

3. 

C, 

0. 

C* 

G. 

0 • 

0* 

0. 

0. 

3* 

c • 

0* 

0. 

G* 

0. 

c . 

0. 

0* 

0. 

3. 

3. 

0* 

Q« 

0« 

Q. 

c • 

0« 

0. 

3* 

3. 

3. 

3* 

c« 

0. 

0. 

0 « 

3* 

•35736r^01 

• 233R5E+n 

0, 

0* 

C. 

0* 

0. 

0. 

D. 

0. 

0* 

Z m 

0. 

G. 

.75035E^D2 

.10505E-*>03 

0. 

c. 

3« 

3* 

3* 

fl 

it • 

3. 

G • 

0* 

0* 

3. 

3. 

3. 

3. 

•46287z^03 

-.11C'^9E«'C3 

0* 

C • 

3. 

0. 

0. 

3 • 

c • . 

3. 

0. 

3. 

0. 

3* 

•2I943E*C2 




0. 

3* 

0 • 

0. 

0« 

3. 

0. 

0. 

0. 

3. 

c • 

0. 

0. 

C • 

c. 

.?• 

0* 

3 « 

c. 

c* 

•12806r^03 

J. 

-•38X61E^03 

-.25503E+C1 

ENTER ETC6P 





0 * 0 « 0 * 0 * 

0 * 0 « 0 * 0 * 

• 35869E+01 .23385E:*01 0. 0. 

0 * 0 « 0 « 0 « 

0 * 0 « 0 * 0 « 

-•29805E:^01 -.22369E:^0 0 


O 



SUHflARY ON THE AR6HENTED HATH IX X 
AR6MENTED MATRIX IS WRITTEN ON LOGICAL UNIT NO. R 
ORDER OF THE COEFFICIENT MATRIX 264 
NO. OF BLOCK ROWS 24 

NO. OF NON-ZERO BLOCKS OF THE ARGMENTEO MATRIX 96 


DECOMPOSE COEFFICIENT MATRIX AND SOLVE LINEAR SYSTEM 


RETURN FROM ETCGP 

AN INDICATION ON THE QUALITY OF THE SOLUTION % 

B IS THE ROW SUM OF THE COEFFICIENT MATRIX «At 

X AND Y ARE RESPECTIVELY THE COMPUTED AND ACTUAL SOLUTIONS OF 
AX=B. THE NORM OF <X-Y) IS •996Q1E-12 


RESMAX= 


406a2£-ll 


IRESM= 24 


JRESM= 12 



HACH NUMBER = .90 


OMEGA = .300 


AW>LITUOE= .026 


****** mode 2f PITCHf PAXIS= 0.000 X/C ****** 

PRESSURE COEFFICIENTS 


I 


BELOW AIRFOIL 


ABOVE AIRFOIL 

ABOVE •• 

BELOW 

X 

9 

C • 

L . 

0. 


.981E+31 

-.383E*01 

-.964E*00 

10 

c. 

3 • 

C. 

^ . 

•790E*01 

-.348E*01 

-.864E*00 

11 

0. 

3- 

j . 

* 

•741E*01 

-*153E*01 

-.723E*00 

12 

a* 

u « 

3. 

0 . 

•597E+01 

•787E*00 

-.526E+00 

13 

3. 

C • 

0. 

r * 

•476£*01 

•14SE*01 

-.250E»BO 

14 

0. 

3 « 

c. 

G. 

•394E*01 

•186E*ai 

•258E-01 

15 

0. 

C . 

c. 

0. 

*617£«>01 

•985£*Q0 

.223E*00 

16 

0. 

G * 

3* 

3. 

.517E*01 

•219E*01 

.364E*00 

17 

0. 

3. 

C. 

0. 

•250E«01 

•326E*01 

.464E*Oa 

18 

0. 

0. 

0. 

3. 

.177E*01 

•310E*01 

.536E*00 

19 

0. 

Q. 

!?. 

?• 

.131E*01 

•289E*01 

.601E*08 

20 

c. 

0 • 

C. 

'■ ^ 

.843E+00 

•250E*01 

.693E«00 

21 

0. 

0. 

C. 

0* 

•AOSE^OO 

.IdSE^Ol 

.821E*0a 

22 

0. 

3. 

c. 

3. 

0. 

D* 

.lOOE^Ol 


TIMEtEtAPSEO TIME <SECONOSI 


1*645 


1.553 



Input for the second sample problem is for a relaxation solution of an airfoil in a Freon flow. 


Card Columns 

No. 1 U 21 31 41 51 ^ 

1 VERSION 1 

2 SAMPLE PROBLEM, AIRFOIL AT M=.9 OSCILLATING IN PITCH. 

3 $PARA FSMACH=0.9, DELTA=0.06, THETA=1.5 , OMEGA=.0600, 

4 GAMMA=1.135, 

7 MSTST=0, 

8 NA=5, 

9 NMAX=200, 

10 ORF=1.8, 

11 NVOL=200, 

12 NP=30, 

13 MODE=2, 

14 PAXIS=0.0, 

15 $END 

Output for the second sample problem starts on the next page. 


43 



vERSIOf.’ 1 


VERSION 1 

RELAXATION SOLUT 10 N *F ULL SPACE 
SU;’SONIC FRECST^U.AM .'NLY 


SAf'PLC PRO'LE^i AIRPOIL AT «t.) OSCILLATIN'. IN PITCH 


79/11/02. «: 


4 ^ 


I 




SPARAH 


FSHACH 

- 


OELTA 

= 

.6E-01. 

THETA 

= 


OMEGA 

= 

•6E-01, 

GAMMA 

= 

.113ST*C1 

AL 

= 

.5£*00t 

IMAX 

= 

25, 

JMAX 

r 

20, 

IS 

= 

0, 

IE 

= 

0. 

JS 

= 

u ^ 

JE 

= 

0, 

NHAX 


200, 

NA 

= 

5, 

ERROR 

= 

.lE-05, 

NP 


30, 

INC 

= 

1, 

ORF 


•18E+21* 

HSTST 

= 

Of 

ISWEEP 

= 

Of 

NVOL 

= 

2G0f 

IP 

= 

Ct 

ILAX 


Of 

URF 

= 

.9E+c:t 

CONPXT 

= 

.12E+C1 « 

CON£:^ 


• lE^O ' « 

MODE 

r 

2f 

PAXIS 

- 

0.-2f 


%ZHO 


MSTSTrO $0 THIS AIRFOIL HAS FIMTt THICK.-^ESS 


THIS IS A ROW R=:lAX OUT-IN 


,-2,S5aC -2.UR2 

-.Ii?58 -*2500 

-1.7383 
. •' ''S3 
1 • 5 - 1 2 

-1.4670 
.2228 
2.0 74C 

-1.2732 
• 3S36 
2.750G 

-1.1347 

.4641 

1.0359 

.5359 

-.9641 

.6012 


-.8636 

.6927 

-.7228 
• 8207 

-&,2';0S -4.3024 

.0620 .1838 

-3. 0656 
.3650 

-2.1177 

.6118 

-1.44C7 

.9571 

-.9571 

1.4407 

-.6118 

2.1177 

— . 365j 
3. 3655 


-.1888 

4.3924 

-.0629 

6.2500 

XK,EPS*XLAM,W1, W2»W7,W9tWll 

.145F*-3i 

.r- 4E*33 

.384E»CC 

.256C*-jO 

.339£<-G0 

.579E- 

Z2 

.2195-31 

.704E+01 

lMAG*y3«W3tyS»W3»W13 c«: 


:.0 .25CF* 

i: 2.0 

.149F-01 O.D 

•256£^20 

Z •'Z 

.3S5E*03 0.0 

.111E*01 


TIMe*£LAPSEO TIME CSrcONnS) 


*>7 4 


TW 




17 


?1 





♦CC 

• 94C f‘ - C 3 

« 4 ♦ y Q 

.357--G2 

.436; *00 

.124F-C1 

• 43S->CC 

.196= -01 

.436E+0C 

• 32.3E-01 

.43'.'-‘*0C 

.357'-01 

.436'*0O 

.402*^-C 1 

.436'+aC 

.419E-01 

.436E-*-uO 

.477r-Gl 

• 4 3 oE *33 

.524*:-Gl 


3 • 


■' * 


TIH£*ELAPSrD TIHi (SECONOS) -O-iS .036 


TIME, ELAPSED TIME (SECONDS) .:81 .CQl 

RELATIVE ERftOR-lST INTEGRAL .0Q03:-» 

RELATIVE . LRR0R-2N0 INTEGRAL .C00337 

RELATIVE EPROR-3RO INTEGRAL .CO::*') 


TIME.ELAPSEO TI>1E (SECONPSJ 



.094 

.013 

TtME.ELAPSEO TI«*£ <SEC0NC.$) 



.097 

.00 5 

TIME.ELAPSEO TIME (SECONnS) 



. *99 

.302 

FARFLO ENTRY 

NMAX = 



TIME, ELAPSED TIME (SECONDS) 



.147 

.048 

N,IERR, JERRfERRMAXl 

I. 

24 

10 

.SRf-Ol 

N,IERR, JERRf IRRM4X1 

2 

22 

11 

.69F-OI 

N,IERR, JERR*-:RRMAX1 

3 

23 

12 

.71^-01 

N,IERR, JERRf^RRMAXl 

4 

23 

13 

.76"-01 

TIME, ELAPSED TIME (SECONDS) 



.275 

.128 

FARFLO ENTRY 

NHAX = 



TIME, ELAPSED TIME (SECONDS) 



.322 

.047 

N,IERR, JERRfERRMAXl 

5 

2 

18 

.94E-C1 

NflERR, JERR tERRMAXl 

f 

24 

IS 

•72E-G1 


.43oE*0C 

.726E-32 


• 269E-: 1 

,43S£*30 

. 583E-: 1 

.43^E*3:. 

*443E- J • 


N * I E!R*^ t JS: KR , r M 4 XI 

T 

24 

1 8 

•8R“-01 

NflERRt J^RR«E:RRMftXl 

R 

24 

1 r 

.81 r-oi 

Nf lERRf JERRiTRRMftXl 

Q 

24 

1-^ 

-75'-31 

TIME, ELAPSED TIME (SECONDS) 


. 4 -^ 4 

.162 

FARFLD 

ENTRY 

NMAX = 

13 


^IHE.ELAP'.'ED TI“E (SFCONOR 

) 


.53: 

.',4 9 

NtlERRi JEIRH*rRRMAXl 

1C 

24 


.71--01 

HiIERRtJERRfERRMAXl 

11 

2 

1 7 

.4ir-Q 1 

NtIERR*JERR*ERRWAXl 

1? 

24 

lY 

.3?s--01 

NtlERR# JERRfERRMAXl 

13 

24 

l4 

.34r.-oi 

NtlERRtJERRtERRMAXl 

14 

23 

7 

.31E-01 

TIME, ELAPSED TIME (SECONDS) 


.8'^2 

.162 

FARFLD 

ENTRY 

N^*AX = 

15 


TIME, ELAPSED TIME (SECOND' 

) 


. 7’ 9 

.04^ 

N,IERR,JERR.ERRMAX1 

Ifi 

24 

- 

.43F-01 

N,IERR, JERRtERRMAXl 

18 

24 

3 

.34F.-01 

N,IERR, JERR.FPRMAXl 

1^ 

24 

4 

.33 -Cl 

N,IERR, JERR,ERRMAX1 

IR 

24 

3 

0 

1 

cc 

N.IERR, JERR,ERRMAX1 

19 

24 

3 

,16r-01 

TIME, ELAPSED TIME (SECONDS) 


.9-9 

• 17C 

FARFLD 

entry 

NMAX = 

2C 


TIME, ELAPSED TIME (SECONDS) 


.^56 

.EAT 

NiIERR«JERRt"RR^AKl 

zc 

24 

•> 

.2TS-01 

NflERRf JERRf ERRMAXl 

?1 

24 

16 

.ISP-Cl 



NtlERRfJEftRt ERR maxi 


22 


2 ^ 


17 


14r-01 


4 ^ 

(D 


N»IERRi JERRfFRRHAXl 23 

JERR*f.RRMAXl 24 

TIHEtELAPSED FIMt ^SECONDS) 

FARFLD ENTRY 

TIMEfELAPSEO TIME tSECONH^:) 


N*IERRtJ£RRf ERR MAXI 25 
N*IERRt JERR^ERRHAXl 26 
N*IERR» JERR.ERRMAXl 27 
N.IERR*JERR*ERRHAX1 28 
N,IERR*UERRfERRMAXl 28 


TIHEtELAPSEO TIM£ (StCONOS) 

FARFLO ENTRY 

TIHEfELAPSEO TIME (SECOND^) 
NflERR.JERRfCRRMAXl 


24 

18 

.12E-ai 

24 

18 

•82E-02 


1.123 .167 


NHAX= 25 

1.169 .046 


24 

19 

•13E-01 

2 

16 

.72E-02 

24 

18 

.68E-02 

24 

7 

.44E-02 

24 

19 

•54E-02 


1.337 .168 

N«AX= 30 

1.383 .346 

24 5 .56E-02 





o 


I 

9FL0W 

AIRFOIL 


A90VE 

AIRFOIL 

8 

-.311F+01 

•149E+C! 


.310E+01 

-.149E+01 

9 

-.314r*Cl 

.15ir-*-Ql 


•313E+C1 

-.151E+01 

10 

-.275F+C1 

•1295+01 


♦275E+01 

-•128E+01 

U 

-.232E+C1 

.138F+n 


•232E+31 

-.13=E*0l 

12 

-,290E+01 

.2125+01 


♦291E+31 

-.2115+01 

13 

-*46Rt>01 

•325F+C1 


•470E+31 

-.321E+01 

14 

-.891'^Cl 

.4GRF+C1 


.9C4E+C1 

-.399E+31 

15 

-.665f +01 

.199r*ci 


•672E+C1 

-.18RE+01 

16 

-.27AE+01 

.26CE+0.; 


•275E+01 

-.176E+0C 

17 

-•2C2E^01 

.49CE-01 


.199E+01 

.136E-01 

18 

-.159<^«-C>1 

-.438F-C1 


•156E+01 

.868E-01 

19 

-.llBE^Gl 

-.IICE^O;? 


.ii4f:+oi 

.123E+0: 

20 

-.737E»GC 

-.139t+Q0 


.70CF+00 

•132E+3: 

21 

-•2QRE+C0 

-.692F-01 


.151E+00 

.41TE-01 

NflERRtJSRRtCRRHAXl 

31 

24 

4 

.44f-02 


NtlERRf JEPRfERRMAXl 

3 2 

24 

3 

•495-C2 


Nf IERRfUERR*t RRKAXl 

33 

?3 

2 

.29E-02 


NiIERR*JERR*5rRR«AXl 

34 

24 

2 

.28F-02 


TIME, ELAPSED TIME (SECONDS) 


1.561 

.173 



FARFLO entry 

NMAX = 

35 


TIMEfELAPSED TIME (SECONDS) 



i.6:'9 

.043 

N.IERRtJERRfERRMAXl 

35 

24 

19 

.32e-C2 

N.IERR* JERR^ERRMAXl 

36 

24 

7 

.175-02 

N,IERR»JERR»ERRMAX1 

37 

24 

6 

.18E-02 

NtlERRf JERRfERRMAXl 

38 

23 

5 

.155-02 

N.IERR.JERR.ERRMAXl 

39 

23 

4 

.135-02 


TIHE, ELAPSED TIME <SECONOS> 1.771 .162 


FARFLD ENTRY 


NMAX= 


A 


AMPLITUDE= .C26 


AHOVE - 

8EL9U 

X 

S21E*:i 

-.298E+01 

-.964E+00 

627E+31 

-.3U2E+C1 

-.8645+00 

550E+31 

-.257E+01 

-.723E+C0 

464E+3 1 

-.273E+C1 

-.526E+C0 

581E»01 

-.424E+C1 

-.2505+00 

958E*-.n 

-.646E+C1 

.2585-01 

130E+32 

-.8G7E+31 

.223E+00 

134E+02 

-.3875+01 

•364E+00 

553E+31 

-.436E+C0 

•464E+00 

4: IE *^01 

-.354E-C1 

.536E+00 

316E+C1 

.131E+00 

•601E+00 

232E+01 

•238E+Q0 

.693E+00 

144E+01 

.271E+00 

•821E+00 

3SOE+50 

•lllE+CC 

.lOOE+01 


041 


t 


TIM£»ELAPSEO TIME (SECONDS) 


NflERR* JERR*ERRMAX1 4C 
Nf lERRf JERR.ERRMftXl 41 
NtlERR* JERR.FRRHAXl 42 
N»IERR» JERRfERRHAXl 43 
N^IERR*JERR.£RR«AX1 44 


TIMEtELAPSED TIME (SECONDS) 

FARFLO ENTRY 

TIMEtELAPSEO TIME (SECONDS) 


N*IERR* JERR*ERRMAX1 4S 
N.IERR* JERRfERRMAXl 46 
N,IE«R»J£RRfERRMAXl 47 
NflERRt JERRtERRHAXl 4R 
N.IERRf JERPtERRHAXl 4^ 


TIMEtELAPSEO TIME (SECONDS) 

FARFLD ENTRY 

TTME»ELAPSED TIME (SECONDS) 


N»IERRf JERR*f RRM4X1 5* 

N,IERR*JERRtERRf^AXl El 

N^IERR^arRR tERRMAXl 52 

N,IERR# JERRfERRMAXl 5T 

N,IERR,JERR.ERRMAX1 54 


TIMEtELAPSEO TIME (SECONDS) 



1.817 

• 04f 

24 

2 

.17"-C2 

o 

5 

.12E-02 

24 

3 

.937-33 

2 

18 

.56:“-03 

24 

15 

.59'-03 


1.970 

.162 

NMAX= 

4^=: 



2*325 

.24', 

2 

.1-5 

.86E-03 

24 

17 

.5ftE-C3 

24 

1ft 

.66E-03 

23 

19 

.40E-33 

24 

1ft 

.32E-03 


2.192 

• 167 

NMAXr 

5: 



2.237 

.347 

24 

? 

.*6f_03 

p 

1ft 

.25"-?3 

■>4 

15 

.26E-.7 3 

23 

16 

.2l"-03 

24 

1 T 

.19t-03 


2.4-:c .163 


55 


FARFLO ENTRY 


N«AX = 


riMCtELAPSCD TIME (SECONDS) 


2.446 


046 


NflCRRt JE^^RfERRMAXl 5*^ 24 19 .21E-03 

N*IERR» JERRiERRMAXl ^6 2 18 .17E-03 

N.IERRt JERR.ERRHAXl 57 24 3 .13E:-03 

Nf lERR.JERRtERRMAXl 5H 2 3 .88E-04 

N,IERRf JERR* "RRMAXl 59 23 18 .83r-C4 

TIH£:#ELAPSED TIME (SECON'^9) 2.6C9 .163 

FARFLO ENTRY NMAX= t,Z 

TIME. ELAPSED TIM-: (SECONDS) 2.655 .346 

NflERRtJERRt* RRHAXl 6C 2 2 .12':'-03 


CJl 

NO 



HACH NUMBER 


OMrQA 





***«*• HOO' ?* PITCH# PAVIS= C#CC0 X/C 


<J^ 


CO 


I BELOW AIRFOIL 


3 

-.311E+01 

.153‘«-01 


9 

-#31?E>G1 

.154-4.T 1 


10 

-♦272^^+:i 



11 

-.229E+C1 

#136E*3l 


12 

-.289F*01 

.214“^01 


13 

-•46EE*C1 

.32RE*:i 


14 

-.8T9r*:i 

.407F+01 


15 

-•654E^ai 

.196E+01 


lb 

-.27ir+ci 

.251E-*«09 


17 

-.196E+C1 

•4C6E-31 


18 

-•154f-»-'!l 

-•49CE-01 


19 

-•112E+01 

-.103E + C-: 


2C 

-.689E+05*' 

-.113E+C0 


21 

-.182f*C" 

-.355r-Ci 


NrlERRtJERR.ERRMAXl 

61 

24 

19 

Nf lERR, JERR^ERRMAXl 

62 

24 

18 

NtlERRt JERRfFRRMAXl 

63 

24 

4 

NflERR* JERRfERRMAXl 

64 

24 

3 

TIHE*£LAPSEO TIME (SECONDS) 


2.833 


FARFLD ENTRY 

NMftX = 

65 

TIME. ELAPSED TIME (SECONDS) 


2.830 

N.IERR.JERR.ERRMAXl 

65 

24 

19 

M.IERR. JERR.ERRHAXl 

Sf 

2 

19 

N.IERR. JERR.ERRMAXI 

67 

24 

6 

Nf ZERR. JERR.ERRMAXl 

68 

23 

5 

N.IERR. JERR.FRRMAXl 

69 

2 

6 

TIME. ELAPSED TIME (SECONDS) 


3.041 


FAPFLD ENTRY 

NMAX- 

73 


060 





ABOVE 

AIRFOIL 

ABOVE - 

BELOW 

X 

.311E+01 

-.153E+C1 

.621E+31 

-.3C5E*:i 

-.964E+03 

.313E^C1 

-.154E*Qi 

.625E^01 

-.308E+01 

-.864E^00 

.272E*:i 

-• 129E+01 

• 544E+-01 

-.259E^01 

-.723E>00 

.229E*C 1 

-.136E*01 

.457E*01 

-.272E*C1 

"-.526E+00 

,289E*c1 

-.214E^01 

.578E^:i 

-.428E^31 

-.250E^CO 

.4S5E*-?1 

-.329E*01 

. 929E^01 

-.657E*- 01 

.258E-01 

.«79E+:i 

-.4G7E*01 

.176E*: 2 

-.S14E+:i 

.223E*>00 

.654E*M 

-.196E+01 

. 131E>C2 

-.392E^C1 

.364E«-0Q 

.2'^2E+01 

-.251E+00 

• 54 3E +0 1 

-.502E*00 

•464E^OO 

.197E^01 

-.4'2"-QI 

. 393E^C1 

-•80SE-C1 

•536E>00 

.154E+01 

•496E-01 

# i^ftE^'Ol 

.985E-01 

•601E«00 

.112E^0l 

.1C6E+Q3 

. 224E+01 

♦211£>00 

.693E+00 

•69CE+0C 

.119£>0‘ 

.138E+01 

.238E«“00 

•821E^00 

.183E>0C 

.36AE-01 

• 364E’*-:0 

.722E-01 

.lOOE^Ol 


•77F-04 


.51F-C4 


17R 


047 

.38E-C4 
• 3S€:-04 
. 26*^-04 

.24F-04 


161 


TIMCiELAPSED T I Mi <r.iCnNO'‘) 




cn 


NflErtR* JCRR 7C 24 *2‘ir-04 
NtlERR » JfRRfERRWiX: 71 2 i .23^-04 
NiIERR*JlRR»lRRMAKI 7? 2 ! .1H"-C4 
N, lERR t JERR tERRMAXl 7-5 2’ .19f-04 
NiIERRiJERRtiRR'^axi *'4 2 2 

TIMEfCLAPStD TIMi (RiCONO'^) 3.2^0 .1^.2 


FftRFLO ENTRY 

riMEfELA^SEO TIHi (SiCON!'?> 


Nf lERRi JERR tPHRMAY! 7*5 

Mt lERR* JERRtERRHAXl 7^, 

Nf lERRf JERRtERRMAXl 7^ 

!V#IERRtJEPRtF RRMAXl 7R 

?4«IERRf JE•■^Rfr:RRMAKl 7i 


TIMEtELA-^SED TIM-: (SECANTS) 


FARFLD ENTRY 

TIMEtELAt>RED TIME (SECON-^^-) 


NtlERP* JERRtERRMAXl R3 

Nf lERRf J*"RR*ERRMAX1 81 

NflERR, JERR tfRRMAXl fir' 

NtlERRf JERRf-:RRM4Xl 8? 

N*IERR* JrRR*‘.RRMA> I 8* 


TIMEfELAPSEO Tn" (SRCONr) 


FARc-lO !"N’’PY 


*J«AX= 

Z.2iS .J4-0 

2 1-' •2C‘-C4 

24 ? .98^-35 

24 5 .9T--:5 

2 .79'>C5 

2 ? .ll"-?4 

3.4 3 .167 

N ■* A Y = 8 J 

5.’.:9 .C46 


7 

1 8 

.ll‘-04 

2 

1 

.9c-;-os 


18 

.Sl'-CS 


1 9 

.47'-:*: 


J 

.41 -3^ 


.1':" 


MMAYr 8’’j 


« 1 


riMEtELAPSEO TIME {SECONDS) 


3*720 


N*IERR» JERRiERRHAKl 85 2 2 

NtlERR* JERRtERRMAXl 3& 2 19 

N,IERR* JERR.iRRHAXl 8? 2 19 

N*IE«R* JERRf^RRMAXl 88 2 17 

l«,IERR, JcRR,ERRHAXl 8P 2 2 

TI.MEf ELAPSED TIME CSrCONOC) 3.885 

FARFLD entry NMAX= 91 

3.933 

2 5 

C71 

CJl 


TIMEtCLAPSED (SECONDS) 

N,IERR» JERRfERRMAXl 90 


4 


• C4S 
•53e-05 
.60E-05 
.46E-05 
.45E-C5 
.31E-05 

• 185 


• 048 
.30E-Q5 



MACH NUMBER = .90 

THFGA = .261 

AHPLITUDE= 

.0 26 

***«.* mode ?♦ PITCH, PAXIS^ 

0.2CC X/C ****** 




o^ESSURE COEFFICIENTS 




I 

H'Loy 

ATRf^OrL 


ABOVE 

AIRFOIL 

ABOVE - 

BELOW 

X 

8 

-.311E*C1 

.153“*-3I 


.311E*:i 

-.157E*01 

.621E+31 

-.335E*C1 

-.964E+00 

9 

-.312'^>51 

.154E+0: 


•312E*01 

•.I54E>01 

.625E+01 

-.ia7E+Gl 

-.864E^00 

10 

-.272E+01 

.12°r+21 


•272E*01 

-♦129E^01 

,544£^01 

-.259£^01 

-.723E^00 

11 

-.229E+01 

.136F+01 


•229r ♦31 

-•136E*01 

. 438£^Q1 

-.272E+21 

-.526E*00 

12 

-.2R9E+01 

.214E+C 1 


.289E*:i 

-.214£*01 

.578E^21 

-.427E*:i 

-•250E+00 

13 

-.465r^G 1 

.328E+01 


.4r>5E^Cl 

-.328F*01 

. 929E+31 

-.657E^0i 

.258E-01 

14 

-.379“+Cl 

•407E*C1 


• 879C*.^1 

-.437E*01 

.176E:^C2 

- .81 3E 1 

•223E^00 

15 

-.654r4-01 

.196-^Cl 


.654F*: 1 

-•196f>Cl 

. 131E+02 

-.392E*01 

•364E>00 

16 

-.272c^01 

.251F*C'' 


.272F+31 

-♦251E+0: 

• 543E-^01 

-.5C1E+DC 

•464E^OC 

17 

-•197E+C1 

.40 lE-C 1 


.197£*C1 

-.4-:0E-01 

. 393E+C1 

-.801E-:i 

•536E+02 

18 

-.154C+01 

- .496E-C 1 


.i:'4£^01 

.496E-:i 

.3C8E^01 

•992E-21 

*601E>00 

19 

-.112F+C1 

-.1C6''-»C-: 


.112E*0 1 

.106E^0: 

.224I4-01 

.211E + :-0 

.693E^00 

20 

-.69:e^cc 



.690r^G3 

.119E^:0 

.158E^01 

.238E+CC 

•821E+00 

21 

-•182r+0Q 

-•364f-Cl 


.182E+3C 

.364E-01 

.3&4E+C0 

.728E-G1 

.lOOE+Cl 

NflERRt JERRf FRRMAXI 

91 

2 

4 

.29^-C5 





NflERRt JERR*£RRMAXI 

92 

2 

19 

.36r-05 





NtlERRf JERRjERRHAXl 

93 

2 

2 

.33^-05 





NrlERRf JERRf ERRMAXl 

94 

2 

19 

•32“-05 





TIME, ELAPSED TIME (SECONDS! 


4.112 

.179 






FARFLO *!' 

NTRY 

NMAXs 

95 


TIME, ELAPSED TIKE (SECONDS) 



4.161 

.C49 

N*IERR» JERRfERRMAXl 

95 

2 

2 

.31F-S5 

Nf lERR* JFRR*£RRMAX1 

96 

2 

1? 

•23E-05 

NflERRt JERR tERRKAXI 

97 

2 

•y 

.25»^-05 

N*IERR* JERRfERRMAXl 

98 

2 

19 

.26^-05 

NtlERR* JERR*ERR«AX1 

99 

2 

y 

.24E-05 


TIMt:*ELAPS£D TIME (SECONDS) 4.^23 .162 


FARFLO ENTRY 


NMAX = 


102 



TIM£*ELAPS £0 TIME (S^.'CONQ'^) 






cn 


N*iERR*JErtR*r*^RMAyi i:" 
N»IERR» JERRfERRMAXl ICl 
Nf lERRf JERRfERRMAXl IC 3 
NtlERR^ JERRtuRRMAXl 1 G 3 
NtlERR* JERRfERRMAXl 104 


TIMEfELAPSEO TIM: {SECONDS) 


FARFLO ENTRY 

TIME, ELAPSED TIME (SECONDS) 


N*IERR*JERR*ERRMAX 1 IC*? 
Nf IERR*JERR#"RRHAX 1 1 C 6 
N*IERR*JERRfERRMAXl 1 C 7 
NtlERRfJERRfLRRMAXI ISR 
N*IERRf JERRtERRWAXl 1 C« 


TIHEfFLAPSEO TIME (SECONDS) 


FARFLO entry 

TIHE#£LAf>SEC TIME (SECONDS) 


Nf rE.RR»J£RR*ERRMAXl U: 

N»IERRtJERR»ERRNAXl 111 
NtlERR^tJERRfERRMAXl 11<7 
NtlERR, JERR*ERR«AXl 113 
N«IERRf JERRfERRPAXl 1 lA 


TIHE*F.LAPSEO TINE (SECONDS) 


FARFLD ENTRY 


2 IV .26F-C5 

2 ? . 24 '- 3 'i 

2 19 •2'E-CS 

2 2 • 19 E- 3 'i 

2 19 .2''*^-25 

.IS"' 

NMAX- ICD 

4.5-13 .C45 

2 ? *1S"-G5 

2 13 .1SE-C5 

2 2 .lSE-35 

2 •17:-C5 

2 2 .15r-D5 

4,7't3 

NNAX= ii: 

4 • T94 .C46 

19 -IS'-Od 

2 ^ . 15 --DT 

2 19 ,14'-15 

? ? .14r-')5 

? 19 .14F-C5 

♦•VbO * 16 ^ 

NNaXr 11 ^ 




TI«E*£:LAPS£0 TI«I (SECONDS) S.CCS 

NtlER'^, JERR^-RRMAXl IIS 2 ? 

M, IERHt JtRR.ERRMAXl 11^ 2 19 

^^,IERR,JERR,LRRH4X1 117 2 2 

N, reRR,jrRR,ERKHAXl 118 2 19 

N^IERRt J£RR*f:RR^AXl 119 2 2 


TIMEtfLAPSEO TIM- (SECONDS) 


5*168 


FARFLD ENTRY NMAXx 12C 

TlHE*"UAPSEO TIME (SECONDS) 5*212 

N«IERR,JERR*ERRMAX1 12? 2 19 


CJl 

00 



045 


• 13E-03 
•12E-35 
.12E-05 

• 12E-05 
•12E-05 

.163 


. 34 * 

• llE-05 





FARFLO ENTRY 

7IHE, ELAPSED TIKE fSECONOS) 


NtlERRiJERR.ERRHAXl 125 
N,IERR,JERR,ERRHAX1 126 
NtlERRf JERR,ERRHAX1 127 
Nf IERR.JERR,ERRNAX1 128 
N.IERR,JERRt£RRMAXl 129 


T1ME*ELAPSED TIRE (SECONDS! 


NRAX= 125 

5.433 .04" 

2 2 .94E-06 

2 19 .91E-06 

2 2 .87E-06 

2 19 .84E-05 

2 2 .81E-0S 

5.593 .168 


FARFLD ENTRY 


NHAX= 13C 


AMPLITUOE= 


026 


ABOVE 

- BELOy 

X 

621£:»01 

-.3G5E*ai 

-.964E+00 

o25E*21 

-.327E^01 

-•864E+00 

544E>01 

-.259E+C1 

-.723E+00 

458E*wl 

-.272E*01 

-.526E^00 

578E4.Q1 

-.427E^0i 

-•250E+00 

929E^0l 

-.657E+01 

.258E-0X 

176E^02 

-.813E^01 

•223E^00 

131E^02 

-.392E*01 

•364E^00 

543E+31 

-.SOIE^OO 

•464£>00 

393E^01 

-.801E-01 

•536E>0D 

30fiE*01 

•992E-01 

•801E>0Q 

224E^01 

.211E«“C0 

•893E^00 

138E*ul 

• 236E«>00 

•821E«00 

364E^00 

•728E-01 

•lOOE^Ol 


TIME, ELAPSED TIME {SE 

CONOS) 


5.640 

• 047 

N.IERR.JERP, ERR MAXI 

130 

2 

19 

.78E-06 

N.IERR, JERRtERRHAXl 

131 

2 

2 

.76E-06 

N.IERR, JERRtERRMAXl 

132 

2 

19 

.74E-06 

N,IERR, JERRrERRMAXl 

133 

2 

2 

.71E-06 

N.IERR, JERR,ERRMAX1 

134 

2 

19 

.68E-06 

TIME, ELAPSED TIME ISc 

ICONOS) 


5-8'3 

.163 

FARELO ENTRY 

NMAX= 

135 


TIME, ELAPSED TIME (SECONDS) 


5.85G 

.047 

N.IERR, JERR.ERRHAXl 

135 

2 

2 

•66E-06 

N.IERR, J’RR.ERRHAXl 

13*3 

2 

19 

.64E-06 

N.IERR, JERR.ERRMAXl 

137 

2 

2 

.62E-06 

N.IERR, JERR.ERRHAXl 

13H 

2 

19 

.60E-06 

N.IERR. JERR.ERRHAXl 

139 

2 

2 

•57E-06 

TIME. ELAPSED TIME (S 

ecoNosj 


6.016 

• 166 

FARELO ENTRY 

NMAX= 

140 


TIME, ELAPSED TIME (SECONDS) 


6.062 

.046 

N.IERR, JERR.ERRHAXl 

14C 

2 

19 

.56E-06 

N.IERR, JERR.ERRHAXl 

141 

2 

2 

.54E-06 

SOLUTION CONVERGEO. 

HAXIHUM ERROR 

IS 

5368E-06 


RESHAX = 

,3C9C6E-Q1 

I = 

24 J = 

2 


TIHE.ELAPSED TIME {SECONDS) 


6.1D3 


091 


MACH NliMPER 





I BELOW AIRFOIL ABOVE AIRFOIL 


8 

-.i82E+a: 

-.364'-ai 

•1R2E*D1 

.364E-31 

9 

-.182F+0: 

-,364"-01 

.lB2c;>02 

♦364E-31 

lu 

-.ispr+o*: 

-.364F-I1 

•ibpe^d: 

.364E-01 

11 

-.182E«-0: 

-.364E-01 

.182E*C? 

.364E-C1 

12 

-.132E + 0I“ 

-.36AE-01 


.364E-01 

13 

-.182"->Cn 

-.364'>01 


•364E-01 

14 

-.182:' + 0: 

-.36Ar-n i 

•182E+0C 

•364E-31 

15 

-.1B2F4-3G 

-.364':-ai 


.364E-31 

16 

-.IBPF^CG 

-.364E-01 

.182F4-0G 

• 364E-n 

17 

-•182r^0Q 

-.364r-Cl 

.182E*2G 

•364E-31 

18 

-.182F-*-GC 

-.364^-01 

.182E-*-aC 

•364E-01 

19 

-.182E+0C 

-.364E-11 

.182F^0C 

.364E-01 

2C 


-.364=--01 

.182£^j: 

•364E-01 

21 

-.ib2f:+gc 

-.364r-0i 

.IB2E^C2 

•364E-01 


TTMEtrLAPSEO TIME <$rCO^iOS> 


6 . 1^5 . 02 ? 



AKPLITUOE= 


26 


A30VE - 

BELOW 

X 

1C3E+C2 

-.314E+01 

-.964E+00 

assE^oi 

-•313E*>01 

-•864E408 

561E^31 

-.257E*01 

-.723E+0C 

462E + 0-1 

-.265E^01 

-.526E^OO 

582E^01 

-.421E*01 

-.250E*0C 

9j6E^D1 

-.638E>01 

•258E-01 

178E+12 

-.82CE4-01 

•223E+0C 

134r>:2 

-.4C0E401 

•364E>00 

545E*C1 

-*438E^GO 

•464E^0G 

394E^:i 

-.143E-C1 

.536E4-00 

3:8E-*'j1 

.167E+D0 

• 601E4-03 

223E*:i 

.280E^0C 

•693E^00 

135E*01 

♦336E+30 

• 821E-»00 


J . 

• lOOE-fOl 
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